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Abstract 

We introduce a unitary dynamics for quantum spins which is an extension of a model 
introduced by Mark Kac to clarify the phenomenon of relaxation to equilibrium. 
When the number of spins gets very large, the magnetization satisfies an autonomous 
equation as function of time with exponentially fast relaxation to the equilibrium 
magnetization as determined by the microcanonical ensemble. This is proven as a law 
of large numbers with respect to a class of initial data. The corresponding Gibbs-von 
Neumann entropy is also computed and its monotonicity in time discussed. 

PACS numbers: 05.30.Ch, 05.40.-a 



1 Relaxation to equilibrium 

The Kac ring model was introduced by Mark Kac to clarify how manifestly irreversible 
behaviour can be obtained from an underlying reversible dynamics, j^JElEl- It explains 
via a simple model some of the conceptual subtleties in the problem of relaxation to 
equilibrium as for example are present in the derivation and the status of the Boltzmann 
equation for dilute gases. In particular, the Kac dynamics shares some basic features 
with a Hamiltonian time-evolution like being deterministic and dynamically reversible. 
In the present paper we extend that dynamics to a unitary evolution on a finite quan- 
tum spin system. Again, the dynamics remains far from realistic but it allows a precise 
formulation and discussion of some features of relaxation to equilibrium for a quantum 
dynamics. That is especially useful and relevant as, in the quantum domain, the prob- 
lem of relaxation is beset with even greater conceptual difficulties. In our framework, 
relaxation to equilibrium becomes visible if one can select a small number of macroscopic 
variables that typically evolve via autonomous deterministic equations to take on values 
that correspond to equilibrium. Typical refers to a law of large numbers with respect 
to the initial data. Paradoxes are avoided by taking serious the fact that relaxation is a 
macroscopic phenomenon, involving a huge amount of degrees of freedom whose evolution 
is monitored over a realistic time-span. One should also keep in mind that relaxation to 
equilibrium goes beyond questions of return to equilibrium, see [2], which are mostly re- 
lated to stability of equilibrium states. A more general introduction to that and various 
related problems can be found in the recent [Sj. 

In sectional we introduce the model and we state the basic result. Section |3 is devoted 
to a discussion of related issues. The proofs are postponed to the final section 0J 
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2 Model and results 
2.1 The model 

Consider TV sites on a ring (periodic boundary conditions). Between any two neighboring 
sites there is a fixed scattering mechanism to be specified below, and at each site, we find 
a spin 1/2 particle. Time is discrete and at each step the ring rotates in a fixed direction 
over one ring segment. Depending on the segment that each spin crosses, it is scattered 
to another state. 

For the Hilbert space Hat we take the TV- fold product of copies of C 2 , the state space at 
each site: 

N 
3=1 

with the standard inner product that defines the Hermitian conjugate for matrices denoted 
by the superscript *. We write the elements oiHN as vectors in (D 2N , denoted by n or rf 
with components r)j indicating the state of the spin at site j. 

2.1.1 Dynamics and observables 

We fix an arbitrary Hermitian matrix H G (D 2x2 , called single site Hamiltonian and we 
construct the unitary matrix 

U = e tH , U* = e~ iH (2.1) 

To each segment on the ring connecting the two nearest neighbors + 1) there is 
associated a scattering variable ej G {0, 1}. The configuration of scatterers is denoted by 
e = (ei, . . . , ejv) and docs not change in time. We define the unitary matrix U € - 3 on (D 2 as 

1 for ej = 

U for ej = 1 

The superscript ej can thus be read as a power and Uj 3 — e lejHj with Hj a copy of H 
working on (D 2 . This defines the scattering mechanism. 

The rotation of the ring is implemented by the operation R on Hn which transforms 
every vector r\— (rji, t]2, ■ ■ ■ , 1}n) into 

Ri] = (777V , ?7i , • ■ • , Vn-i) G U N 

and R* = R^ 1 , rotation in the opposite direction. Combining that with the scattering 
mechanism finally gives rise to the unitary operator Un = Un{c) on Hn via 

JV 

U N = R(g) 

3=1 

To be specific we look at positive times and the state at time t = 1, 2, ... is then obtained 
from the state rj G Hn at time t = from 




rj t = U N (t) r? = U N (t - l)U N tj 
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This unitary dynamics on the finite-dimensional TLn is rather simple. For many observ- 
ables (like the total magnetization) the dynamics is entirely equivalent with leaving the 
spins in place and rotating instead the scatterers. We therefore work with the more 
convenient 

N 

U N {t) = ^uf (t ' t] (2.2) 

with kj(e, t) = 53 n =i e j-n (modulo TV). The updating of the spins is independent modulo 
the fact that they may have a scattering mechanism in common. 

The Hamiltonian is a sum of single site contributions, 

N 

H N = Y,Hi (2.3) 

with Hj = 1®...®!?®...®! (on the j-th site) and for concreteness, we decompose 
the single-site Hamiltonian in the Pauli basis 

H = hi a x + h 2 <r y + h 3 u z = h ■ a ^ h G R 3 (2.4) 

The eigenvalues arc e + = \\h\\ = h and e_ = — e+. (We can of course ignore adding 
a constant to this Hamiltonian.) Our observables are also built up by adding one-site 
contributions. We are given a Hermitian matrix A S (D 2x2 and we construct the system 
operators 

JV 

A W =^1®...®A®...®1 (2.5) 

where the A is at the j-th position. We will consider the magnetization vector Mn = 
(Mjy, Mff, M^) defined through the one-site observables M a — a a , the Pauli matrices, 
with a € {x, y, z}. 

The restriction of looking only at one-site observables will be discussed in Section l3~4l 
2.1.2 Equilibrium 

The dynamics depends non-trivially on the scatterers but the equilibrium properties are 
independent of them. That is similar to the situation in the Boltzmann-Grad limit where 
the hard core matters dynamically but does not enter in the computations of energy or 
pressure. 

Energy is conserved in the sense that the Hamiltonian commutes with the time evolution: 
[U N (t),H N ] = 0. 

Equilibrium is characterized by the microcanonical distribution. Let V'+iV'- be the two 
eigenvectors of H with eigenvalues e + > e_ and spectral projectors P e+ and P e _ respec- 
tively. For arbitrary e £ [e_, e+] we select an energy space by 

e e 

where, in we sum over all (ej) with fixed Y2j e j satisfying Ne < J2j e j < e + — e_ + 7Ve; 
the normalization Z„ ensures that Tr[pN] = 1. Alternatively, we could have summed in 
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(|2.6|) over all (e^) with ^ . ej/N in a certain interval around e and at the very end let the 
interval shrink to zero. The average 

{A N }? = Tr[P? A N ] (2.7) 

for a Hermitian matrix An on TLn defines the (finite volume) microcanonical ensemble. 

The (infinite volume) equilibrium magnetization corresponding to an energy e is then 
defined from 1)2.7)1 to be the limit 

1 ~ N 

"^Jfe 1 iWe (2-8) 

This limit exists and can easily be computed: it equals 

m e = e A (2.9) 

2.1.3 Initial data 

Initial data are determined by a density matrix p^ on and by the choice of the 
scatterers e. Concerning the state of the spins, let us keep in mind an initial preparation 
with a particular magnetization in the z-direction, not in equilibrium. Let Q+ and Q- 
be the projectors on spin up and spin down respectively, with <j z = Q+ — Q— , to define 

Q^ = J2Q mi ®-..®Qm N me{-l,+l} (2.10) 

m 

where, as in 1)2.6)1 . in Y] m , we sum over all (rrij) with fixed ^ . rrij satisfying Nm ^ 
J2 3 rrij <2 + Nm. Equation (|2~TH| gives the projector on the magnetization space with 
magnetization converging to m as N f +oo. A possible initial density matrix is then 
(0 < d(m, N) < +oo is a normalization) 

"»" = «^o« < 2 '"> 

Another (but thermodynamically equivalent) choice would be the grand-canonical 

Po = exp (\M* N ) (2.12) 

and similarly for the other directions of the magnetization. A more general class of initial 
density matrices p$ will be introduced at the beginning of section Most important 
is that they satisfy a law of large numbers for observables of the form l|2.5|l and that the 
single site marginal remains well-defined in the thermodynamic limit, see l|2.22|l . 

Now the dynamics starts and the density matrix at time t is defined from (|2.2(l : 

p^=^(e) = U N (t)/^U N (tr (2.13) 
and now also depends on the scatterers ey. 

We are interested in the magnetization in p^ 1 and in its limiting behavior as (first) ./V f +oo 
and (then) t f +oo. In order to take these limits, we also need to specify the initial (and 
unchanging) condition of the scatterers. For this we introduce the set 



lee{0,l} N 



N 1 

J2tj = \»N]\ (2.14) 

3=1 ) 
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with \pN] the smallest integer not smaller than fj,N. The constraint fixes the fraction 
of active scatterers to be about p. The probability to select one particular e in f2j^ is 
uniform: 

1 



(2.15) 



All limits will involve that probability distribution. We say that a sequence of functions 
Gat on typically takes the value g^, written 



e- lim G N = 



iff 



lim M«[(G N - 9li )v]=0, p=l,2 



(2.16) 



where the expectation refers to the probabilities (|2.15|l . Via the Chebyshev inequality one 
can reformulate (|2.16|) as a statement about the fraction of scatterers in for which 
Gn — g^- The last ~ is made stronger if higher p's in 1(2.16(1 are obtained. For simplicity, 
we restrict us here to the average (p = 1) and to the variance (p = 2). 

2.1.4 Entropy 

One expects for the given model that the one-particle distribution satisfies an autonomous 
equation. That information is encoded in the single site density matrix. For finite N, the 
marginals are defined through 



,N 



(2.17) 



with Tr 3 the reduced density matrix at site j. The Gibbs-von Neumann entropy is defined 
following the ideology of [3 through the variational principle 



S N (p N )=su P -Tr[p> N \ogp> N ] 

p ,N 



(2.18) 



where the supremum is over all density matrices p lN with marginal v N from l|2.17|l . 
Obviously, the supremum in l|2.18|l is attained for the product state p lN = (^) v N and 
therefore through 1)2.17(1 . the quantum entropy equals 



S N {p 



-NTt 



v N logv N 



We insert from (|2.13() in 12.19fl and put, hoping all goes well, 

A/J+oo JV 



(2.19) 



(2.20) 



as the time-dependent entropy (density). The more frequently considered von Neumann 
entropy -Tr [p^ logp^] does not change in time and is therefore here less relevant. 



2.2 Results 

The conditions for our main result on the dynamics ((2.1 3|) are on the level of initial data, 
see section 12.1.31 We ask that p$ is a density- matrix that is dispersionfree in the sense 
that for all observables 1(2. 5|) and for all continuous / 



lim 

AfT+oo 



Tr 



-/ Tr 



An n 



= 



(2.21) 
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Furthermore, we need the existence of the thermodynamic limit of the system-averaged 
one-site marginal defined in (|2.17l) : 



1 N 



1 + mo • a 



j'=i 



and similarly, we also define the time-evolved version of vq, i-e., 

1 N 

Vt = e- lim — >^ Tr, fori 



Finally, the sequence Pq must satisfy a technical condition. We write 

r(N) N 
j— 1 



(2.22) 



(2.23) 



(2.24) 



the initial density matrix as a sum of products of 2 x 2 matrices and we assume that 

x (2.25) 



sup sup \\xf,w\\= C < 

7Vg]N ™ = l,..,r(N) 
j = l ; ..,N 



These three conditions are obviously satisfied by l|2.12Jl . That (|2.21|l is also satisfied for 
(|2.11|) needs an argument given in the example after the proof of Lemma 14.31 

The initial magnetization is denoted by mo = Tr[oVo] and the equilibrium magnetization 
rh e is defined in (|2.8|l . We denote components of the magnetization as m a , a 6 {x, y, z}. 
The initial energy e is found from e = Tr[i/Vo]. 

Theorem 2.1 

For all continuous functions f and for initial Pq satisfying the above conditions, 



e- lim Tr 

N]oo 



[Ml 

V N 



Pi 



= /«) 



with 



rh t = m e + Re 



toq — me + 



i(h x Toq) 



1 — /i + fj, exp[i(e + — e_) 



In particular for f{x) = x, the magnetization 



lim e-lim — Tr \Mn p?~\ = m e 

tT+ooWT+oo N 1 J 



(2.26) 



(2.27) 



(2.28) 



converges exponentially fast to its equilibrium value, for all \i G (0, 1) when e+ — e_ is not 
a multiple of 2tt. 



The entropy H2.20J1 equals 



l + \\mt\\ , 
«(*) = o lo S 



1 + 



l-ll^l 



log 



1- ||TO t | 



(2.29) 



where m t is given by (|2.27|) . s(t) is increasing in time (H— theorem) and can be written 
as 

s(t) = S[v t ] , S[u]=- Tr[is log v] (2.30) 
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3 Discussion 

3.1 Kac ring model 

The Kac ring model in 6 , p. 99, can be recovered by the appropriate choice of the Hamil- 
tonian (|2.4|l . by only looking at the magnetization in the z-direction. Identifying rji = 1 
or |t) with a white ball and the state rji = —1 or \\) with a black ball in the original ring 
model, the switching between a white and a black ball corresponds to a spin flip. This is 
accomplished by the single-site Hamiltonian (written out in the a z base): 

which, up to an irrelevant constant, corresponds to the choice h = — (■§, 0, 0) in l|2.4|l . The 
eigenvalues of this Hamiltonian satisfy e+ — e_ — n and the equilibrium magnetization is 
m% = 0. We start from randomly sampling the balls (or spins) with a fixed overall color 
(or magnetization) in H2.11[l so that mf — nv( = for all t. Substitution in H2.27[l yields: 

m z = m z (l-2p t ) t 

That coincides with the relaxation formula in the original Kac model. It is interesting to 
observe that the relaxation to equilibrium gets slower for all the extensions that we have 
considered here. 

3.2 Molecular chaos 

The Stofizahlansatz, or repeated randomization hypothesis consists in replacing the real 
dynamics by an effective dynamics in which the scattering mechanism is not kept fixed 
but gets replaced with an average. In that way, memory is being erased of where the 
scatterers are. An effective dynamics then works on the single particle level and since 
the fraction of active scatterers equals p, the one-site density matrix v t is either copied 
with probability (1 — p) or is replaced with a scattered density matrix with probability 
p. More precisely, 

T*(v) = (l-fi)i/ + fiUuU* (3.1) 

defines a completely positive map (the quantum equivalent of a discrete time stochastic 
dynamics, see £Q). The unitary U should be substituted from i|2.1|) . On the level of the 
full density matrix p N , we just make the product of copies of T* over all N sites so that 
a product state is mapped into a product state: If v N (j) = Tij [p N ] , then 

N 

A*(p N ) = <g)T*(v N (j)) (3.2) 

That is the dual of a map A and we write A t for the map A applied t times. The effective 
dynamics thus amounts to replacing the time evolution U^it) ■ Ujsr(t)* (depending on the 
scatterers €j) by A£ = (^T^ (only depending on p) and expectations at time t become 

Tr[A t (A N ) pg] 

Computations are here even simpler. 
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The results of section 12.21 say essentially that the above effective dissipative dynamics 
reproduces the correct result (see in fact Lemma f4. 211 . For example, 

e-lim ^Tr [Mj^f] = km I Tr [A t (M£) p?] = Tr [r t (a> ] (3.3) 

In other words, the unitary dynamics l|2.2|) (depending on the scatterers ej) is typically 
equivalent with the effective dynamics T t on the one-particle level. The marginal v t of 
(|2.22|) exists and can be obtained from v t +\ — T*v t . That specifies the equation (|2.27(l 
and determines H2.30|l . 

Note finally that the reduction of the unitary time evolution B2.2JI to an effective dynamics 
on a one-particle system remains non-Abelian. Of course, our macroscopic observables 
start to commute, 

[ N ' N 1 \N 

but r* remains defined on density matrices and does not yield a classical dynamics on 
diagonal elements of the density matrix Pq. 

As stated in Theorem 12. II as in the original Kac ring model (or as for the rigorous deriva- 
tion of the Boltzmann equation), one really does not need this assumption of molecular 
chaos to get relaxation. It is replaced with statistical assumptions on the level of the ini- 
tial conditions for the scatterers — they are typical with respect to the probability l|2.15(l 
- and for the spins — they are e.g. randomly sampled with a fixed magnetization in 
(|2.11|) . In that case we think of both the r\ (wavef unction) and the e (classical scattering 
centers) as dynamical variables. 

3.3 Autonomous equations and H— theorem. 

The equation (|2.27|) is autonomous in the sense that the value of rht determines fht+i 
once we know h = hri = (e + — e_) n/2 in the Hamiltonian l|2.4|l and the fraction // of 
scatterers: 

m t +i = frit — 2/t [(H x rht) sin h cos h — h x (n x rh t ) sin 2 h] (3-4) 

A more suggestive expression is obtained by decomposing rh into a parallel and perpen- 
dicular component along h: 

rh = (rh ■ n)n + h x (rh x ft) = rh^ + 

for which 

m\\(t + 1) = m\\(t), fh±(t + 1) = (1 — 2/isin 2 h)rh±(t) + 2/isin/icos/i rhj_(t) x n (3.5) 

The relaxation formula l|2.27|l simplifies for certain initial conditions. For example, when 
taking H2.11fl or l|2.12|l . so that mg = ttiq = 0, we get the damped oscillator 

rrif — nig + (toq — m^r* cos ujt 

with r 2 = (1 — /i) 2 + /i 2 + 2/i(l — /i) cos 2h and tan(w) = [i sin 2h/ (1 — fi + p cos 2h). The 
frequency is maximal for h = tt/A and the damping is maximal for h — 7r/2. In that case, 
the relaxation of ml also looks autonomous but that is because of the special initial data. 
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In general, does not relax autonomously and in contrast with l|2.29[) . its associated 
entropy 

1 + mf , 1 + ml 1 — mf , 1 — mf 
log log - 

2 5 2 2 5 2 

need not be monotone (even in the case of (|2.11Jl or (|2.12(l 1. Going to the bigger picture 
with three magnetization components reveals a new structure in which the entropy does 
increase. From (|3.5l) .we see that ||m(i)|| decreases as 

\\m ± {t + 1)|| 2 = (1 - 4/i(l - fi) sin 2 h) \\m ± (t)\\ 2 (3.6) 

and hence, combined with (|2.29(l . the monotonicity s(t + 1) > s(t) follows and strictly 
so if fj, S (0, 1) and h ^ kit. We thus see here that also in the quantum case, autonomy 
on some macroscopic scale is intrinsically connected to monotonicity of the corresponding 
entropy 

3.4 Relaxation to equilibrium? 

The relaxation can be read off clearly from (|3.5|l : we get a spiral motion in the plane 
perpendicular to ft. The model does however not show the full glory of relaxation to 
equilibrium. That is already true in the original Kac ring model. As an example, consider 
a macroscopic variable which involves a two-spin function, in contrast with the single site 
observables that we introduced in 1)2.5(1 : 

3 

It is easy to verify that A^/N does not show relaxation. This is not surprising given 
the fact that two neighboring spins live exactly the same history and thus they do not 
decorrelate. This seriously restricts the usefulness of the model for studies of relaxation 
but it remains possible to study the phenomenon of relaxation for the special macroscopic 
observables of (|2.5(l . 

On the other hand, the fact that we consider relaxation in terms of macroscopic observ- 
ables should not be considered as a restriction. After all, equilibrium is characterized by 
a maximal entropy condition given macroscopic constraints. In particular, relaxation is 
not be read off from the Liouville-von Neumann evolution for all microscopic details of 
the density matrix. The fact that we first let N f +oo (before time) avoids the presence 
of (quasi-)periodicities or of Poincare-recurrence. 

A final useful comparison concerns the notion of ergodicity which is quite universal, see 
e -g- .5, 4. Consider a quantum mechanical system with a discrete non-degenerate energy 
spectrum (E n ) with (4> n ) a complete set of orthonormal energy eigenfunctions. The wave- 
function at time t is denoted by ip(x,t) = exp[— iHt]ip(x) for initial ip( x ) = ^2 n c n4'n{x). 
The time-average of the expectation value of a Hermitian A is 

A = lim ^ / dt { dxip*{x,t)Aip(x,t) 
TT+oo i Jo J 

It is rather easy to see in that case (but it remains true in a much broader context) that 



A = ^2\c n \ 2 / dxcj)* n (x)A(j) n (x) 

n 



(3.7) 
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which could be argued to correspond to the microcanonical ensemble. Note however that 
no use was made in the above of the fact that N is large nor of the fact that we consider 
macroscopic variables. The result l|3.7|l is equally not very satisfactory: we are interested 
in the typical manifest behavior over realistic time-spans for large systems while (|3.7|) 
gives trivial information about an infinite time-average that cannot be identified with 
relaxation phenomena. There is in fact no relaxation to equilibrium on a microscopic 
scale even though l|3.7|) always holds. 

3.5 Higher spins 

Instead of considering spin 1/2 particles, it is also possible to set up exactly the same 
problem as above for a Hilbertspace TLn — ®^Li ■ We are then dealing with n x n ma- 
trices but there is little difference in the computations except for one extra complication: 
there is a larger class of conserved quantities. To explain that, we split the single-site ob- 
servable A (a Hermitian matrix in <C nxn ) in a part that commutes with the Hamiltonian 
H, and its orthogonal complement: 

A = Ah ® A± (3.8) 
with [Ah, H] = and with respect to the scalar product 

(A\B) = Tr [A*B] (3.9) 
The set of commuting observables Ah is a vector space with dimension 

dim(„4#) = n 

If n > 2 there is more than just the 1 and H G Ah (e.g. not all functions of H can be 
expressed as linear combinations of 1 and H). 

Because of the triviality of the dynamics as was discussed in the previous section, these 
extra conserved quantities give rise to conserved quantities for the full system dynamics 
and none of them disappear as the size of the system grows. For more realistic systems, 
truly interacting, we expect that, in the thermodynamic limit, essentially only such phys- 
ical quantities as particle number and total energy remain conserved. In that respect, 
the case n = 2 is more physical and that is the reason why we have concentrated on it 
since the beginning. The mathematical theory can however be completed if we change the 
definition of the microcanonical ensemble (|2.7Jl to include the extra conserved quantities. 

4 Proof of results 

Since it is important for physical interpretation that the relaxation curve i|2.27[l is obtained 
for almost all choices of scatterers e drawn from 12.14JI . we must first remind in Lemma 
14.11 of the thermodynamic equivalence with a grand-canonical set-up. Next, in Lemma 
14.21 is shown that an effective dynamics reproduces the time-evolution of our macroscopic 
observables. Using that, and from Lemma 14.31 the propagation of condition 12.21fl is 
obtained. It is therefore sufficient to derive 12.27fl to get l|2.26[) . That computation is 
introduced by Lemma 14.41 The entropy 1|2.29J) and its monotonicity can then also be 
computed using the obtained law of large numbers. 
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Let fl = {0,1} and e 6 Cl . We have already denned the "canonical" measure P^; in 
(PTT5|) on of ipnijl . Define now its "grand-canonical" version, the probability mea- 
sure P M on f2 N , as the Bernoulli measure with density P AI (e J = 1) = fi. E^(-) and E M (-) 
denote the expectation values with respect to P^ and P^ respectively. The following 
equivalence is standard. 

Lemma 4.1 For all functions g = <?(ei, ■ ■ . , £&) on fl k 

2 2k+1 k„ 



E»-E M ( ff ) 



iV 



i/ 2fc < N 



with 



Proof: 



\g\\ = max \g(£)\ 



E^-E^) 



< E iff(oi 1^(0-^(01 



The right-hand side can be bounded by developing P^ (£) in conditional 
expectations. By explicit calculation 



1^(^+1 = 1 1 (6,..-, &) = 0)-H < 

after which a simple calculation finishes the proof. 



fe + 1 
N-k 



We define for fixed j and matrix element per element 

r t (^) = E„^(i,e) = e-lim E*A,-(t,e) 

iv J+oo 



where the last equality follows from the previous Lemma l4.ll 



Lemma 4.2 For macroscopic observables An = Ylj Aj and continuous functions f, the 
dynamics is typically equivalent to the semigroup dynamics A t ,i.e., 



e-lim Tr 

JVT + oo 



N > Pt 



= lim Tr 



/A, 



A N 
N 



Proof: 



The proof is for all monomials. By the Stone- Weierstrass theorem, this proves the 

theorem for all continuous functions /. 

Let 

D(n) = {l,...,N} n 



and 



T(n,t) = {if G D[n)\ Vi,j < n :i ^ j => \K i -K j \>t\ 
A{K) = A)§ ® . 



.4 



k(0 

m 
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where 

K = {j(l),...,j(l),...,j(l),...,j(l)} 
fe(i) fc(i) 
We also need the index set J(if ) = {j \ Bp : K p = j'}. 

We begin by calculating the e-average. Denote Ajy(t, e) = UN(t)*A^UN(t), and 
similarly for the one-site observables Aj. 



Tr 



TV 



w KGT(n,t) jeJ(K) 
1 

TV™ 



KeD(n)\T(n,t) 



n] 



The last term can be dropped because 



^ E Tr[^)(M)^] <^ E WMK)\\=0 

K£D(n)\T(n,t) K e D {n)\T (n ,t) 



u\\ n 

N 



We apply Lemma PTTI with gK,w{t) = Iljeif ^ [A?'(*' e )Xjfm] as function on fi tn . 
For if e T(n,t), 



< II Tr [Mt,e)X?, w ] - (E^Tr [A,-(t,e) 



We replace E M Tr[Aj(i, e)xjfj by Tr[r t (A,)x"J and we add terms of order iV" 1 
(essentially the ones we first subtracted) to arrive at the average from the lemma. 

For the variance, we similarly write 



A N 
N 



Po 



N 2r, 



E K ( Tr W0(*> e X] Tr W L X*> e)Po] ) + 

if,Z,GT(ra,t) 

Also the contribution of sets K and L that are within a distance less than t are of 
lower order. We apply Lemma 14. II hist as above. 



Lemma 4.3 For all continuous functions f and macroscopic observables A 



N- 



e-lim Tr 

JVT + oo 



A T 



N \ M 



N 



PI 



= lim f Tr 

AfT+oo \ 



A, 



— \o N 



at all times t. 
Proof: 



By Lemma l4.2l we can insert the initial density matrix. The statement then follows 
from applying condition H2.21[l to the observable ■ T t (Aj). 
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As an illustration of the lemma above, we show that it can be applied for the state Pq 
given by H2.11JI . 

Example: 

Call P^ the projector on the eigenvalue a of the macroscopic variable Am/N . Let 
7 be the spectrum of Am/N and I = UI N . We want to establish a large 
deviations property for the probability measures W N : 



W N {K) = 



Tr [P a N 



a£l N nK 



(4.1) 



The product P^ Pq can be written as a sum over 1-dimensional projections and 
(|4.1(l can be developed and written in terms of 

Tr[P a .Q m .] 

which takes four possible values depending on (pj, nij). The problem amounts 
therefore to estimate the probability that at the same time ^2jPj, Y]j vn,j and 
Y]j rrijPj take on specific values when the (j>j , rrij ) are independent and identically 
distributed with uniform weights. That can be found from the multinomial 
distribution and Stirling's formula gives the large deviation rate function 
corresponding to l|4.1|l . As a consequence, a law of large numbers is satisfied and 
in the thermodynamic limit the state l|2.11[l concentrates on only one eigenvalue of 
A N /N. 



We now come back to (|2.22|l defining the marginal v$. We split the Pauli matrices a a — 
a% + cr" as in (|3.8|1 and parameterize in the single-site Hamiltonian basis: 



p p 
1-p 



o e a 
e a o 



The typical value of the magnetization is defined to be: 



m t = e- lim Tr 

TVT+oo 



M N n 

IT* 



(4.2) 



Lemma 4.4 

rht = m e + Re 



/ _ _, i(h x mo) 
I mo - m e H 



Proof: 



(l -/i + /iexp[i(e + -e_)]) 



Fixing i, we write k{e) = Ylj=i e j ana ^ we obtain 



(4.3) 



Z a {t) 

z a (t) 
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with 

z Q (t) = 6 a (l - /i + /icxp[i(e+ - e_)])' (4.4) 
<jfl is invariant under the time-evolution and it is easy to verify that 

lim Tr [r t (a Q )i/ ] = Tr [a%v ] = m a e 

This yields 



2f39 a {l — (/, + /iexp[z(e + - e_)]) 



mf = m° + Re 
Now it is a matter of writing 2/39 a as 



2f30 a = Tr[is a° A 

ii 

then inserting H — h ■ a, and recognizing the resulting complex numbers as 
components of the vector product h x mo, to find (|4.H|> almost immediately. 



The asymptotics (|2.28(l is trivial and the proof of Theorem 12.11 is thus completed if we 
add the equation l|2.29|l for the entropy S[v t ]. 



Proof: 



The explicit formula (|2.29|l can immediately be obtained by using the Bloch 
representation: 

1 + m t • tr 

** = 2 

Inserting it in the entropy density l|2.30|l yields the result. 
Concerning the increase of entropy: denote by v e the one-site marginal of the 
equilibrium state H2.6(l . with magnetization rh e . A short computation gives 

p 
1-p 

where p is the first diagonal element of v$ , see 114. 2|) . 
The relative entropy between two marginals is defined as 

S[v\v e ] - TrHog^] -Tr[^log^ e ] (4.5) 

We immediately have: 

S[vt\v e ] = S[v e ] - S[v t ] 

since the T-dynamics works only on the off-diagonal elements of vq. Together with 
the entropy-contraction inequality for completely positive maps 

S[P> t )|P> e )] ^S[vt\v e ] 

and since T*(i> e ) — v e , we find: 

s[r k (vt)]=s[v t +i]>s[v t ] 
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